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ABSTRACT 

We report a bimodality in the azimuthal angle distribution of gas around galaxies as traced by Mg II absorp- 
tion: Halo gas prefers to exist near the projected galaxy major and minor axes. The bimodality is demon- 
strated by computing the mean azimuthal angle probability distribution function using 88 spectroscopically 
confirmed Mgll absorption-selected galaxies [Wr(2796) > 0.1 A] and 35 spectroscopically confirmed non- 
absorbing galaxies [W,-(2796) < 0.1 A] imaged with HST and SDSS. The azimuthal angle distribution for 
non-absorbers is flat, indicating no azimuthal preference for gas characterized by Wr(2796) < 0.1 A. We find 
that blue star-forming galaxies clearly drive the bimodality while red passive galaxies may exhibit an excess 
along their major axis. These results are consistent with galaxy evolution scenarios where star-forming galax- 
ies accrete new gas, forming new stars and producing winds, while red galaxies exist passively due to reduced 
gas reservoirs. We further compute an azimuthal angle dependent Mgll absorption covering fraction, which 
is enhanced by as much as 20-30% along the major and minor axes. The Wr(2796) distribution for gas along 
the major axis is likely skewed toward weaker Mg II absorption than for gas along the projected minor axis. 
These combined results are highly suggestive that the bimodality is driven by gas accreted along the galaxy 
major axis and outflowing along the galaxy minor axis. Adopting these assumptions, we find that the opening 
angle of outflows and inflows to be 100° and 40°, respectively. We find the probability of detecting outflows is 
~60%, implying that winds are more commonly observed. 

Subject headings: galaxies: halos — galaxies: intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

It is well estabhshed that Mgll AA2796,2803 ab- 
sorption, detected in background quasar/galaxy spec- 
tra, arises from gas associated with foreground galax- 
ies and thus provides a unique means to directly ob- 
serve mechanisms by which galaxies acquire, chemically 
enrich, expel, and recycle their ga seous component (see 
IChurchill. kacprzak. & Steidell l2005l for a review). The 
Mgll ion is ideal since it traces metal-enriched low-ionization 
gas wifli 10'^ < Af(Hi) < 10^^ cm" ^ (IChurchill et al.ll2()00t 
iRigbv. Chariton. & Churchilll 120021) and, as a result, is 
detected out to projected galactic ra dii of 100 kpc 
(iKacprzak et alJ2008tTchen et al.ll20 10aV A significant quan- 
tity of Hi gas is probed by Mgll absorption in galaxy ha- 
los with roughly 15% of the gas resid ing in DLAs and 5% 
of the total hyd rogen in stars (.Kacprzak & Churchilll 120 iTt 
iMenard & Fukug ita 2012). 

A large body of evidence suggests that Mgll ab- 
sorption traces both outflo ws from star- fo rming 
galaxies ("Bouche et al.' "200^, 'Tremonti et al.' "2007; 
Zibetti et a l. 2007; Martin & Bouche 2009; Weine ret al . 
2009; Chelouche & Bowen 20 101; iNestor et all boflT 
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2010a 



gal ax- 
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Ribaudo et al 20 11; . Kacprzak et al. 2012; Rubin et al. 2012; 
Martin et al.! 120121) . WTiile it is clear both processes are 
occurring, it is difficult to disentangle which absorption sys- 
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tems may be uniquely associated with either process. Since 
outflows are presumably metal-enriched, whereas accreting 
material should be metal-poor, comparing the absorption-line 
and host galaxy metallicity provides one suita ble test for 
discriminating winds from accretion (e.g., iRibaudo et al.l 
201 1; Kacprzak et al. 2012, and references therein). However 
it is not yet feasible to perform this experiment for a large 
sample of galaxies. 

The Mg II spatial distribution relative to their host galaxies 
also provides a promising test. Outfl ows are expected to ex- 
tend along the galaxy minor axis (e.g.. lStrickland et al.ll2( )04). 
whereas accretion would progres s along fi laments more 
planer to the galaxy (e.g., Stewart et al,|l201_l'). Modeli ng the 
morphologies of 40 HST imaged absorbers, Kacprza k et al.l 
(,2011b) found that the rest-frame Mgll A2796 equivalent 
width, W, (2796), is dependent on galaxy inclination, suggest- 
ing a co-planer geometry. Furthermore, the Mg II kinematics 
are consiste nt with be ing couple d to the gala xy angular mo- 
mentum dSteidel et aLll2002l; IKacprzak et al.ll2([Hdal) . These 
results provide evidence that some Mg II absorbing gas is ac- 
creting co-planer to galaxies. 

Stacking ^ 4000 background galaxy spectra. lBordoloi et al.l 
(|2()11) showed that, statistically, unresolved Mgll-doublet 
equivalent widths are larger along the galaxy projected minor 
axis than along the projected major axis within ~ 40-50 kpc, 
suggesting that the strongest absorption is ejected within 
winds along the minor axis. Using 10 absorbers. lBouche et al.l 
( 120121) showed that Mg II absorption arises along the projected 
minor and major axes of the host galaxies, suggesting that 
Mgll is primarily detected in outflows and accretion, how- 
ever, the inclusion of a control sample of non-absorbers is 
necessary to vali date the bimodal claim . Both works were fur- 
ther validated bv lChurchill et al.l (l2012h who, using 65 galax- 
ies, found that Mgll absorption arises preferentially along 
the galaxy major axis and then along the galaxy minor axis 
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as Wr(2796) increases; they also found that for Wr(2796) < 
0.1 A, the galaxy inclinations and azimuthal angles with re- 
spect to the quasar sight-line are consistent with random dis- 
tributions. 

In view of these results, we aim to further explore the az- 
imuthal distribution of Mg II absorption for a large sample of 
spectroscopically confirmed galaxies. In this Letter, we com- 
pute the mean azimuthal angle probability distribution func- 
tion for galaxies hosting Wr(2796) > 0.1 A absorption and 
show that it is bimodal, peaking near the galaxy minor and 
major axes. We also show that the distribution is flat for galax- 
ies with Wr(2796) < 0.1 A (i.e., non-absorbers). We further 
study this distribution in terms of galaxy colors and equivalent 
widths. We adopt a /i = 0.70, JIm = 0.3, 57a = 0.7 cosmology. 

2. THE SAMPLE AND ANALYSIS 

We compiled a sample of 88 spectroscopically confirmed 
Mgll absorptio n-selected gala xies with W,(2796) > 0.1 A 
[33 from Kacprzak et al.l (12011b.) imaged with HST, 46 from 
iChen et all (I2010aft and 9 from iKacprzak et al.l (1201 lal) im- 
aged with SDSS] and 35 spectroscopically confirmed non- 
absorbers with W,(2796) < 0.1 A [21 from iChurchill et al.1 
(l20T2h imaged with HST, 14 from lChenetal.1 (12010 a') im- 
aged with SDSS] having 0.09 < z < 1.1 ((z) = 0.33) and 
9 <D < 200 kpc ((D) = 48 kpc). The galaxies were selected 
to be isolated; group and double galaxies are not included. 
The Wr(2796) = 0. 1 A cut maximizes the number of absorbers, 
while retaining a significant number of non-absorbers. 

Galaxy AB magni tudes and colors were determined follow- 
ing the methods of iNielsen. Churchill. & Kacprzak! (120121) . 
Rest-frame B-Kor B-R were determined from the measured 
apparent magnitudes. When only one color was measured, we 
applied a conversion determined from a linear least-squares fit 
to the average c olors from ga l axy sp ectral energy distribution 
templates of M annucci et alJ (120011) . We obtained rest-frame 
colors for all but one galaxy. 

For this work, all galax y morphological properties were ob- 
tained using GIM2D dSimard et alj|2002l) . 56 galaxies were 
imaged with HSTPNFPC2, for whic h we adopt the rn orpholo- 
gies and orientations measured by Kacprzak et al. (zOlT9) 
and IChurchill et al. (2012). For the Kacprzak et al. (2011a|) 
galaxies imaged with SDSS, we a dopt the GIM2p morph olo- 
gies and orientations measured bv lKacprzak et al.l (1201 lah and 
iBouche etal] (120121) . 

The SDSS galaxies taken from ' Chen et al.l (l2010ah have a 
median redshift of (z) = 0.25. Given the short SDSS image 
exposure times and the high redshifts, we studied the reliabil- 
ity of the GIM2D models. We identified 13 SDSS galaxies 
that had also been imaged with HST from our previous work 
((z) = 0.50). All SDSS galaxies in our sample are resolved 
with >20 pixels at 1 .5 cr above the background in the r-band. 
Using the SDSS r-band images, we computed azimuthal an- 
gles using GIM2D and then also from the 1.5 a isophotes 
using Source Extractor (Berli n & Arnoutsll 19961) . In order to 
minimize the number of GIM2D free parameters, we fitted the 
data using an exponential disk only, with zero bulge contribu- 
tion. 

We found that both the GIM2D and isophotal azimuthal 
angles of the SDSS galaxies are consistent with the GIM2D 
azimuthal angles of the HST galaxies (within 1 a uncertain- 
ties of the HST values). In general, the errors on the GIM2D 
SDSS values are a factor of a few larger, whereas it is diffi- 
cult to quantify uncertainties for the isophotal values. Thus, 



we modeled the SDSS r-band imaged galaxies as an exponen- 
tial disk with GIM2D, applying the modeling prescription of 
IKacprzak etal] (1201 lal) . 

We modeled all galaxy morphological types since 86% of 
early-types have regular stellar rotation (14% exhibit slow ro- 
tation), spanning the full range of apparent ellipticities, with 
90% being kinematically and photometric ally aligned (see 
lEmsellem etani201 UlKrainovid et all201 Ih . 

We adopt the convention of the azimuthal angle <i> = 0° to 
be along the galaxy major axis and <I> = 90° to be along the 
galaxy minor axis. 

Direct binning of the azimuthal angles, the approach taken 
bv lBouche et al.l f2012l ). applies only when the measured un- 
certainties are smaller than the bin size. Further complica- 
tions arise when the uncertainties are asymmetric, as is the 
case for GIM2D model parameters. Here we take a different 
approach: we represent the measured azimuthal angles and 
their uncertainties as univa riate asymmetric Gaussians (see 
iKato. Omachi. & Asd 120021) . thus creating an azimuthal an- 
gle probability distribution function (PDF) for each galaxjQ. 
From the continuous azimuthal PDFs, we then compute the 
mean PDF as a function of The mean PDF represents 
the probability of detecting Mgll absorption at a given 
This technique provides higher weight per azimuthal angle 
bin for galaxies with well determined $. However, even the 
less robustly modeled galaxies provide useful information; 
the method is equivalent to stacking low signal-to-noise spec- 
tra or images to search for a coherent signal. 

3. THE BIMODAL AZIMUTHAL DISTRIBUTION 

In Figure [T] we present the binned mean azimuthal angle 
PDF for the 88 absorbing and 35 non-absorbing galaxies. The 
binned PDFs have been normalized to the total number of 
galaxies in each sub-sample. The shaded regions about each 
bin are the 1 a deviations produced by a jackknife analysis in 
which we randomly removed 10% of the sample for a million 
permutations. 

The non-absorbers exhibit a relatively flat distribution, as 
expected for random <I>, and the fluctuations are consis- 
tent with Poisson no ise. This result confirms the work of 
IChurchill et al.l (120121) . who showed that the azimuthal angle 
distribution is random for W, (2796) < 0.1 A. 

The absorbing galaxies exhibit a bimodal <i> distribution, 
with one peak corresponding to the background quasars prob- 
ing the galaxy projected major axis and a second, higher and 
broader peak, when the galaxy is probed along the projected 
minor axis. The shape of the mean PDF indicates that Mgll 
absorbing gas is preferentially located along the projected ma- 
jor and minor axes of galaxies. Our result confirms the inno- 
vatory work of Bouche et al. (2012) based upon 10 galaxies. 
A plausible scenario is that wind/outflow gas is distributed 
about the galaxy minor axis with relatively high frequency, 
whereas accreting/infalling gas is frequently found in a co- 
planer geometry. For the remainder of the discussion, we re- 
fer to major axis absorption as infalling or accreting gas and 
to the minor axis absorption as outflowing or wind gas. 

We further divided our sample of absorbing galaxies into 
two impact parameter sub-samples, those with D > 40 kpc 
and D < 40 kpc, and found that the bimodal shape of the mean 
PDF is well preserved for both impact parameter sub-samples . 

PDFs using the quadratic model of IBarlowl 420031) produce similai' re- 
sults, however their PDFs contain large unrealistic spikes when the uncer- 
tainties are small. 
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Fig. 1. — (top) Binned azimuthal angle mean probability distribution func- 
tion for 88 absorbing galaxies (solid line) and 35 non-absorbing galaxies 
(dashed line). The areas of the histograms are normaUzed to the total number 
of galaxies in each sub-sample. The shaded regions are 1 a confidence in- 
tervals based upon a jackknife analysis in which 10% of the sub-sample was 
removed at random. Given galaxies exist in multiple bins, the size of the Icr 
confidence intervals does not simply scale as the square-root of the number of 
galaxies per bin. However, the sizes of the confidence intervals will increase 
if more galaxies are removed in the jackknife analysis. The non-absorbers are 
consistent with a random distribution within the expected Poisson noise. The 
bimodal distribution suggests a preference for Mg II absorbing gas toward the 
galaxy major axis ($ = 0°) and along the minor axis ($ = 90°). — (bottom) 
Azimuthal dependence of the covering fraction with shaded regions provid- 
ing the Icr confidence intervals. Note the higher covering fractions along the 
major and minor axes. 

This implies that outflowing and accreting gas maintain their 
major and minor axis orientation preferences, respectively, far 
out into the halo. 

The width of each peak in the mean PDF provides rough 
constraints on the geometry of outflowing and inflowing Mg II 
gas. The peak at <I> = 0° suggests that accreting gas is found 
within A<I> ~ ±20° of the galaxy plane. For $ = 90°, the 
peak is initially narrow, suggesting the majority of galaxies 
expel gas within A$ ~ ±10° of their minor axis. The much 
broader tail of the $ = 90° peak suggests that the opening an- 
gle of outflowing gas may be as large as A$ ~ ±50°, but also 
the frequency at which this gas is detected decreases for larger 
opening angles. Our estimated outflow opening angle is con- 
sistent with other published values determined using different 
techniques (e.g.. .Walter et al...2002: .Martin et al...201Z) . The 



15 



10 - 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



I 



J L 



_l I I I I I I I I I L_ 



20 40 60 

$. Azmuthal Angle 



80 



Fig. 2. — Galaxy color dependence of the azimuthal distribution of Mgll 
absorbers. The solid line (bl ack) is the d i stributi on shown in Figure [T] The 
color selection adopted from IChen et aTl j2010al) of S — i? < 1.1 represents 
late-type galaxies (dashed blue line) and B — R> 1.1 represents early-type 
galaxies (dotted red line). Shaded regions are Icr confidence intervals based 
upon a jackknife analysis. The data are consistent with star-forming galaxies 
being dominated by outflows. 

ratio of the area between the outflowing and infalling gas, bi- 
furcated at 40°, suggests that ~ 60% of Mgll absorbing gas 
is outflowing. 

In Figure [T] we plot the Mgll gas covering fraction, which 
is the ratio of the number of absorbers to the sum of absorbers 
and non-absorbers in each azimuthal bin. This first presenta- 
tion of the covering fraction as a function of azimuthal angle 
shows a peak at 0.8-0.9 along the projected minor axis, a de- 
crease of 20-30% at intermediate $, and an increase again 
toward the projected major axis. The mean covering fraction 
of our sample, 72%, is consistent with previous results (e.g., 
iKacprzak et alilOW^ IChen et alMKoS . 

3.1. Colors, Equivalent Widths, and $ 

We separated the sample of absorbers into early-type and 
late-type g alaxies based on their rest-frame B-R colors. Fol- 
lowing Ch en et al.l (l2010al) . we adopted the color cut B-R < 
1 . 1 to represent late-type (blue) galaxies and B-R > 1.1 to 
represent early-type (red) galaxies. In Figure|2] we present the 
mean azimuthal angle PDF broken down by galaxy color The 
early-type galaxies exhibit a relatively flat distribution with 
~ ±1 (T fluctuations across bins. However, the elevated fre- 
quency at $ = 0° is suggestive that gas is accreted along the 
major axis, which may provide the coo l gas reservoirs ob - 
served around early-type galaxies (e.g., iGrossi et al.ll2009b . 
The late-type galaxies clearly dominate the overall bimodal 
distribution. This result is intuitive if we speculate that red 
galaxies have significantly less accreting gas, and therefore 
little fuel for star formation, resulting in no outflowing gas. 
For blue galaxies, the accretion is higher, providing fuel for 
star formation that then produces outflows. We also split the 
absorber sample by galaxy color to yield equal numbers of 
red and blue galaxies in each sub-sample; the same trend ap- 
pears. The azimuthal bimodality of blue galaxies and the flat 
distribution of red galaxies is consistent with the picture in 
which blue star-forming galaxies exhibit outflows due to ac- 
creted gas fueling star formation. 

One might argue that red galaxies may have a preference 
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Fig. 3. — (left) The area-normalized distribution of ^,-(2796) for blue galaxies with D <40 kpc, blue galaxies with D> 40 kpc, and red galaxies . The weighted 
mean of each distribution is shown with the solid points above the histogram; the eiTors are equivalent to the point size. While consistent with'Bordoloi et ajj 
1201 II) . we found blue galaxies at lower D having larger W,-(2796) compared to red galaxies and blue galaxies at D > 40 kpc. — (right) The W,-(2796) dependence 
for 3> > 40° and $ < 40° separated into two D bins. The W,-{219(s) weighted mean for each sub-sample is shown above the distributions with $ > 40° indicated 
by circles and 4> < 40° by arrows. The data suggest a difference in the W,-(2796) distributions as a function of both D and <&. This may imply that either there is 
less gas column density along the projected galaxy major axis then along the projected minor axis, or, the gas could be metal-enriched along the galaxy minor 
axis compared to the major axis. This would be expected in a accretion/wind scenario. 

sub-samples. Similarly, for $ < 40°, the mean W,.(2796) dif- 
fer by a factor of 3.8 for the smaller and larger D sub-samples. 
The Wf(2796) distribution for outflowing gas at smaller D 
is characterized by a paucity of Wf(2796) < 0.5 A absorp- 
tion and a higher abundance of Wr(2796) > 1 .5 A absorption, 
whereas the W, (2796) distribution for outflowing gas at larger 
D is characterized by an abundance of Wr(2796) < 0.5 A ab- 
sorption and a paucity of Wr(2796) > 1 .5 A absorption. A KS 
test yields that the populations differ at the 2.6a significance 
level. 

Comparing the D < 40 kpc W,.(2796) distributions for wind 
versus accreting gas, we note a similar trend; at smaller im- 
pact parameters, accreting gas exhibits a higher abundance of 
WX2796) < 0.5 A absorption and a paucity of WX2796) > 
1 .5 A absorption. KS tests yield that the outflowing gas with 
$ > 40° and D < 40 kpc differs from the inflowing gas with 
$ < 40° and D > 40 kpc at the 3. la significance level. Quali- 
tatively, these combined results may imply that either there is 
less column density and/or velocity spread of gas in projection 
along the galaxy major axis than along the minor axis, and/or, 
the gas may be more metal-enriched along the galaxy minor 
axis than along the galaxy major axis. These trends would be 
expected for the accretion/wind scenario we have described. 



for w eaker absorption (see IZibetti et al.ll2007t iBordoloi et alj 
120111) . We thus investigated whether flie Wr(2796) dishibu- 
tion depends on galaxy color and D. In Figure [3] we show 
the area-normalized distribution of W, (2796) for red galax- 
ies and for blue galaxies with D < 40 kpc and D > 40 kpc 
(we have too few red galaxies to separate them into im- 
pact parameter bins). The blue galaxies at smaller D have 
a higher optimally weighted mean Wr(2796) than blue galax- 
ies at larger D. This is consistent with the well-known anti- 
correlation b etween Mgll equivalent width and impact pa- 
rameter (cf.. iNielsen. Churchill. & Kacprzak|l2012 , and ref- 
erences therein). The sub-sample of red galaxies has a mean 
W,.(2796) consistent with the blue galaxies at larger D. 

Interestingly, Figure|3]shows a paucity of Wr(2796) < 0.5 A 
absorption at smaller D for blue galaxies compared to both 
blue galaxies at larger D and to all red galaxies. Also, note the 
paucity of Wr(2796) > 1 .5 A absorption at larger D for blue 
galaxies as compared to both blue galaxies at smaller D and to 
all red galaxies. However, a KS test yields that the Wr(2796) 
distribution of blue galaxies at smaller D are statistically con- 
sistent with the Wf(2796) distributions of both the red galaxies 
and the blue galaxies at larger D (2.0 a, a high significance 
level is difficult to achieve with small number statistics). 

The suggested weaker absorption associated with red galax- 
ies and with blue galaxies at larger D in our sa mple is consis- 
tent with the findings of Bordoloi et al.l (1201 ll ). who interpret 
this as winds extending out to projected distances of ~ 40- 
50 kpc. However, the bimodality in the azimuthal distribution 
is present for our D > 40 kpc absorber sample, suggesting 
winds persist beyond 40 kpc, but with smaller Wr(2796) indi- 
cating that the wind gas thins out at larger D. 

We further investigated whether the distribution of 
W,(2796) differs for infalling gas or outflowing gas by split- 
ting the D < 40 kpc and D > 40 kpc absorbers into the two az- 
imuthal bins, $ < 40° for inflowing gas and $ > 40° for out- 
flowing gas. The value $ = 40° is the point of inflection in the 
mean azimuthal angle PDF. In Figure [3] we present the area- 
normalized Wr(2796) distributions. For $ > 40°, the mean 
Wr(2796) differ by a factor of 3.3 for the smaller and larger D 



4. CONCLUSION 

We have demonstrated that the distribution of Mg II absorp- 
tion around galaxies exhibits an azimuthal angle bimodality 
when compared to non-absorbing galaxies, whereby Mg II ab- 
sorption is preferred near the projected major and minor axes. 
Outflows and accretion likely extend beyond 40 kpc, but with 
smaller mean Wr(2796). Furthermore, the blue star-forming 
galaxies drive the bimodality, suggesting that the accretion of 
gas drives star-formation that produce outflows. Red galaxies 
may exhibit some gas accretion along the major axis. Under 
these assumptions, we compute opening angles of outflows 
and inflows to be 100° and 40°, respectively. 

In addition, the data suggest that star-forming galaxies at 
low impact parameters exhibit higher mean (2796) than red 
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galaxies and star-forming galaxies with larger impact param- 
eters. Also, the Wr(2796) distribution for the infalling gas 
along the major axis exhibits a smaller mean Wr(2796) than 
outflowing gas along the minor axis. This could result from 
less gas being probed along the major axis than the minor 
axis, less velocity structures/dispersions expected for winds, 
or higher metal enrichment along the minor axis, as expected 
for wind driven material. The probability of detecting out- 
flows is '^60%, implying that winds are more commonly ob- 
served, likely because the opening angle of outflows is 2.5 
times larger than for accreting gas. 

The data paint a picture that is consistent with the general 
idea of how the buildup and evolution of galaxies occurs as 
well as how galaxies enrich their gaseous halos. Our results 
signify that additional HST imaging of absorbing and non- 
absorbing galaxies will substantially increase our understand- 



ing of galactic scale feedback and accretion of intergalactic 
gas. 
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